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Lasercooled RaF as a promising candidate to measure molecular parity violation 

T.A. Isaev, 1 ^ S. Hoekstra, 2 -0 and R. Berger 1 '^ 

1 Clemens- Schopf Institute, TU Darmstadt, Petersenstr. 22, 64287 Darmstadt, Germany 
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The parameter W & , which characterizes nuclear spin-dependent parity violation effects within the 
effective molecular spin-rotational Hamiltonian, was computed for the electronic ground state of ra- 
dium fluoride (RaF) and found to be one of the largest absolute values predicted so far. These calcu- 
lations were performed with the complex generalised Hartree-Fock method within a two-component 
(quasi-relativistic) zeroth-order regular approximation framework. Peculiarities of the molecular 
electronic structure of RaF lead to highly diagonal Franck-Condon matrices between vibrational 
states of the electronic ground and first excited states, which renders the molecule in principle suit- 
able for direct laser cooling. As a trapped gas of cold molecules offers a superior coherence time, 
RaF can be considered a promising candidate for high-precision spectroscopic experiments aimed 
at the search of molecular parity-violation effects. 
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Diatomic radicals containing heavy nuclei are, together 
with closed shell chiral molecules with high Z nuclei [l|,|2|, 
especially suitable for the search of parity violation effects 
which depend on the nuclear spin (see review [3(). Up to 
the present time, the only measured nuclear-spin depen- 
dent parity violating (NSDPV) quantity is the anapole 
moment [J, |5| of the 133 Cs nucleus J6J. Nevertheless, 
it was concluded that the result of this experiment is 
in contradiction with previously obtained restrictions on 
NSDPV constants (see review [7j), which calls for further 
attempts to determine NSDPV effects. Currently, an ex- 
periment on measurement of the nuclear anapole moment 
of barium in the barium fluoride (BaF) molecule using a 
molecular beam is under preparation at Yale University 
[8| . Here we propose the use of trapped radium fluoride 
(RaF) molecules as a promising candidate to probe parity 
violation, which has never been detected in molecules. 

In recent years a number of techniques has been devel- 
oped to prepare, control and cool molecules [9j. A level 
is now being reached where these techniques can have a 
significant impact on precision experiments. Stark decel- 
eration, well suited to slow down and trap diatomic polar 
molecules, has been developed in the last decade |10l - tl3 | 
and the possibility of laser cooling of molecules has only 
very recently received more attention |13l4l5| . If such 
cooling methods can be successfully applied to sensitive 
molecules, a new generation of precision measurements 
can be done. 

As we show in this paper, the RaF molecule appears 
very well suited for laser cooling methods, because of 1) 
the highly diagonal Frank- Condon matrix between the 
electronic ground state and the lowest electronically ex- 
cited state as well as 2) the facts i) that the correspond- 
ing electronic transition frequency is in the visible region 
with a reasonable lifetime and ii) that isotopes with nu- 
clear spin quantum number 1—1/2 are available. Com- 
bined with its large W a , this makes the RaF molecu- 
lar a prime candidate for a precision experiment probing 
molecular parity violation. 

An experiment with radioactive molecules requires a 



specialized facility. At the KVI of the University of 
Groningen the TRI/iP facility on the production of ra- 
dioactive isotopes was commissioned in the year 2007, 
with the perspective of using these radioactive isotopes 
in precision experiments to test fundamental symme- 
tries [l6(. As the accuracy in a spectroscopic mea- 
surement is ultimately limited by the coherence time, 
there is a clear potential for precision experiments us- 
ing trapped species. A number of precision experiments 
using trapped species has been set up, such as Ra + to 
probe parity violation [ljj and Ra atoms for electron 
electric dipole moment (eEDM) detection [18|. We will 
discuss the production, deceleration and trapping of RaF 
molecules furtheron in this paper. 

Calculation of the W a parameter Previously [191 ] we 
utilised the quasi-relativistic two-component (2c) zero or- 
der regular approximation (ZORA) for calculating the 
NSDPV parameter W a of the effective molecular spin- 
rotational Hamiltonian in a number of heavy-atom di- 
atomic molecules. The corresponding P-odd contribution 
to this Hamiltonian reads as W a A;>i[AxS' c ff] •/ with the ef- 
fective electron spin operator S b q, the nuclear spin opera- 
tor /, the vector A which points from the heavy nucleus to 
the light one and the parameter Ha of the heavy nucleus 
7] that includes the effect of the nuclear anapole mo- 
ment and the NSD Z° exchange contribution. One may 
additionally attempt to absorb effects from the nuclear- 
spin independent Z° exchange contribution in fc^, but 
as this contribution depends via the hyperfine coupling 
tensor also on the specific electronic structure, it may ap- 
pear more canonical to include this contribution in the 
electronic structure parameter W a . For the time being, 
however, we have omitted this (supposedly small) addi- 
tional term. 

The results obtained with this approach are at the level 
of paired generalised Hartree-Fock (paired GHF, for the 
classification used see review |20() for off- diagonal matrix 
elements (MEs) between the degenerate electronic states 
while diagonal MEs correspond to complex GHF, thereby 
including spin-polarisation effects. The advantage of our 



approach is that it can be straight-forwardly extended to 
polyatomic molecules. In the case of diatomic molecules 
with 2 S ground state, we can obtain reliable estimates 
of the NSDPV effect together with the hypcrfinc cou- 
pling tensor A usin g ZO RA. In our calculations we used 
a modified version [H [H HI of the TURBOMOLE 
code 23[ [24 J with an even-tempered atomic basis set on 
Ra, a basis set of triple-zeta quality from the TURBO- 
MOLE basis set library augmented with diffuse p- and 
d-functions on F atom (see supplementary material at 
[URL will be inserted by AIP] for basis sets and other 
parameters of the calculations) . A finite spherical Gaus- 
sian nucleus model is used and a ZORA model potential 
as prop osed by van Wiillen |25j with additional damp- 
ing [26j was applied. The equilibrium distance R e of the 
electronic ground state was obtained from complex GHF- 
ZORA calculations and found to be in good agreement 
with R c determined in later four-component (relativis- 
tic) Fock space coupled cluster singles and doubles (FS- 
CCSD) calculations. As one can see from the Table U 
the changes for all parameters are within 6 % when the 
R c value is taken either from FS-CCSD or GHF-ZORA 
calculations. The computed absolute value of the W a pa- 
rameter is in the kHz range and thus one of the largest 
predicted for diatomic molecules so far. As Ra nuclei 
have an even number of protons, measurements of W a 
in RaF would complement ideally previous 133 Cs exper- 
iments on the nuclear anapole moment of nuclei with 
an odd number of protons. The hyperfinc coupling ten- 
sor calculated herein for RaF with accounting for spin- 
polarisation effects is highly isotropic and can be used, 
once measured, for estimating the quality of our predic- 
tion for W a . 



TABLE I: Molecular parameters for the 2c-GHF-ZORA de- 
scription of the electronic ground state 2 £i/2 of 223 RaF and 
225 RaF. The results are given for both R c taken from the 
GHF-ZORA (4.40 a ) and FS-CCSD (4.24 a ) calculations. 
The parallel and perpendicular components of the hyperfine 
tensor are given in GHz, Wa in kHz (with estimated depen- 
dence on the Ra mass number being on the order of a few 
Hz). 

223 Ra p 225 Ra p 

-R c /ao | A||/GHz Ax/GHz I Ay /GHz Ax/GHz | HykHz 



4.40 



4.24 



1.90 



1.82 



1.86 



-15.1 



-15.4 



-14.8 



-15.1 



1.36 



1.30 



Besides the high W a value there is also another pecu- 
liarity which makes RaF very attractive for experimental 
NSDPV search, namely that this molecule appears well 
suited for being directly cooled with lasers. 



Franck- Condon factors and direct laser cooling of 
molecules One of the major problems in direct laser 
cooling of molecules is the spontaneous radiative decay 
into a manifold of rovibrational states not covered by 
the cooling cycle. Transitions between rotational states 
are driven by selection rules for vector operators and can 
be reduced to the two- level case [l4j, but depopulation 
of vibrational levels during cooling cycles is the obsta- 
cle. In Ref. 13[ several molecules were proposed which 



have highly diagonal overlap matrices between vibra- 
tional eigenfunctions of the ground and the first excited 
electronic states (in terms of the molecular spectroscopy 
— highly diagonal Franck-Condon matrices). For 10 4 to 
10 5 cooling cycles, the probability of the molecule to stay 
inside the cooling loop has to be greater than 0.9999 with 
the fourth digit being significant. In practice this means 
that one has to know FC factors with an accuracy of 
0.01%, which is extremely challenging to achieve in mod- 
ern calculations (although FC factors can in principle be 
measured with high accuracy). Thus predictions of FC 
factors with the accuracy required by experimental needs 
appear practically impossible. 

Nevertheless, one can try to narrow down the num- 
ber of candidate systems for direct cooling with laser 
based on peculiarities of the molecular electronic struc- 
ture which renders certain molecules promising. We have 
analysed the electronic structure of diatomic molecules 
with one electron over closed shells. Such molecules have 
large magnetic moments (in comparison with the closed 
shell systems) which makes them suitable for magnetic 
trapping. More detailed considerations shall be made in 
a later article, here we would just like to give the short 
summary. Highly diagonal FC matrices should be ex- 
pected when the unpaired electron is making transitions 
between orbitals that do not contribute to the bonding. 
Such one-electron wavefunctions can be of the following 
type: 1) "Electron in lone orbital". This situation is 
well-known in theoretical chemistry, when, e.g. as a re- 
sult of sp™ hybridisation, an unpaired electron is directed 
out of the bonding region. This takes place e.g. in the 
electronic ground state of MF and MH, where M belongs 
to the group of alkaline earth metals (Mg to Ra). Two 
molecules from this series were already used to indicate 
the perspectives of direct molecular cooling, namely CaH 

2) "Atom-like electron 
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in Ref. [13[ and SrF in Ref. 
in molecule". In this case, the outer electron is, due to 
some reasons (symmetry as a rule), located in a highly- 
symmetrical state, which has practically atomic nature. 
Such a situation happens for example in HI + [27| , where 
the unpaired electron resides in an orbital which is mainly 
of atomic P3/2 character centered on iodine. In isoelec- 
tronic neutral molecules (such as TeH) the situation is 
expected to be the same. 3) In heavy-atom compounds, 
the valence electron orbital is often quite diffuse, thus 
the maximum density of the unpaired electron is located 
further away from the molecular core than in lighter ho- 



mologues. Thus among all possible molecules, those with 
heavier nuclei appear promising. 

The heavy-atom molecules RaF, HgH, HgF etc. which 
should belong to the first class ("Electron in lone or- 
bital") and HI+, PoH, TeH etc. which should belong to 
the second class ("Atom- like electron in molecule") are 
expected to have highly diagonal FC matrices and, thus, 
appear, in principle, highly suitable candidates for direct 
cooling of molecules. Unfortunately, the molecules be- 
longing to the second class are not particularly suitable 
for NSDPV experiments as the unpaired valence electron 
is in an orbital with high total angular momentum and 
does not have noticeable spin density in the vicinity of 
heavy nucleus. As to the first class, the molecules HgH 
and HgF have (as RaF) a 2 E ground state and, accord- 
ing to our estimates, even higher W a values (for HgH it 
is 3.3 kHz, although we would like to emphasize again 
that spin-polarisation and other electron correlation ef- 
fects are not accounted in our W a estimates, while these 
effects can well be important). The excitation energy of 
the first excited state in these molecules, however, is in 
violet or ultraviolet region (for HgH 407 inn. for HgF 25G 
nm) [28J. Especially at 256 nm it is difficult to produce 
sufficiently intense laser light. In the series MF, instead, 
where M belongs to the group of alkaline earth metals, 
the excitation energies from the ground to first excited 
electronic states are in the visible region, so one can ex- 
pect a similar situation in RaF. This wavelength region 
is directly accessible using diode lasers. 

To substantiate the considerations above, we per- 
formed exemplary four-component Fock space (FS) 
CCSD calculations of the six lowest-lying states of RaF 
using the DIRAC08 package [2J|. The details of the FS- 
CCSD approach are given in the review [30|. Electrons 
of the 6s,6p,7s shells of Ra and 2s, 2p of F are correlated, 
the basis sets and other computational parameters can 
be located in the supplementary material. 

The two energetically lowest lying electronic states are 
identified as 2 Ei/2 (ground state) and 2 H\/2 (first excited 
state). Then the points calculated in the interval 3.5 ao to 
7.0 ao were fitted by a Morse potential (the points and 
fitting curves are given in Figure [1]) and spectroscopic 
parameters together with the corresponding FC factors 
were determined. The FC factors were calculated by the 
Mathematica script for FC factors calculation between 
vibrational Morse oscillator states taken from Ref 
The results obtained are presented in Table HH 

The laser cooling transition ( 2 n 1 / 2 (w'=0) <— 
2 Si/ 2 (w=0)) is found to lie around 710 nm, and an 
upper state lifetime on the order of ~ 25 ns is expected. 
With these parameters the capture range for optical 
molasses is ~ 5 m/s, and the doppler limit temperature 
is ~ 150/j.K. If a precooled sample of molecules can 
be created, as is described in the next section, only 
a few thousand photons would have to be scattered 
per molecule to reach this temperature. At these 
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temperatures only a very small volume in the trap is 
occupied, opening the possibility to perform precision 
measurements using the trapped molecules. 

The sum of the FC factors for the ground vibrational 
state of the 2 Ii-\/2 state and three vibrational states of 
2 Ei/2 state is equal to 0.9999 with four digits after a 
point being stable in respect to variations of the fitting 
parameters of Morse potential (see supplementary mate- 
rials for details). We would like to emphasize here, that 
high digital precision in calculations of FC factors does 
not mean high accuracy of the spectroscopic constants 
calculations. Reliable estimate of the sum of FC factors 
which can be used for practical purposes is expected to 
have two digits after a decimal point. Nevertheless both 
numerical results and arguments from electronic struc- 
ture theory provide, in our opinion, good reasons for con- 
sidering RaF for measurements of NSDPV effect and for 
direct laser cooling. 



TABLE II: Estimated molecular spectroscopic parameters for 
the electronic ground state 2 £i/ 2 and the first electronically 
excited 2 iI 1 / 2 state from FS-CCSD calculations of RaF. The 
dissociation energy D was obtained as the difference between 
total energies of the states at the equilibrium distance R e and 
at 9.0 ao without attempting to correct for basis set superpo- 
sition errors. Fitting parameters of Morse potential with the 
corresponding standard errors can be found in supplementary 
materials. Parameter uj c is the standard harmonic vibrational 
wavenumber. We would like to emphasize that the accuracy 
of the determination of the parameters is lower (especially for 
D and T c ) than the number of digits reported. Nevertheless 
we report the results with high digital precision to provide 
reaper points for future calculations. 



R c /ao We/cm 1 D c /cm 1 To/cm 



2 E 1/2 4.24 428 32104 



2 n 1/2 4.24 432 31302 14012 



Proposed measurement strategy and sensitivity to NS- 
DPV signal The W& parameter can be obtained in a di- 
atomic molecule by a measurement of the parity-violating 
coupling between hyperfine-rotational levels of opposite 
parity. An experimental procedure to do so in a molec- 
ular beam has been proposed in Ref !8| , where tuning of 
the close-lying spin-rotational levels of opposite parity to 
near-degeneracy due to the Zeeman effect is used to in- 
crease the P-odd asymmetry in the measured signal. The 
accuracy in the determination of W a that can be reached 
in such an experiment is proportional to the count rate 
of the molecular signal (dN/dt) , the interaction (coher- 
ence) time (T) and the total measurement time (r) as 
SW a - l/(Ty/TdN/dt). 

For an interaction region of 5 cm and a molecular beam 



velocity of 500 m/s, the interaction time is on the order of 
100 /its. As it has been demonstrated that cold molecules 
can be trapped for at least a number of seconds 32|, [33| , 
the potential gain in accuracy is 4 orders of magnitude. It 
is also clear that part of this gain is lost as the number of 
molecules will probably be low; so production and deliv- 
ery of the RaF molecules to the working volume is crucial 
for the proposed measurements. On the other hand the 
high Franck-Condon overlap allows for the detection of 
multiple photons even from a single trapped molecule. 

The two radium isotopes that are most suitable are 
213 Ra and 225 Ra, with lifetimes of 164 s and 14.9 days, 
respectively. Both isotopes have a nuclear spin I = 1/2, 
which reduces the number of lasers required for lasercool- 
ing. The 213 Ra isotope has been successfully produced 
and used in experiments at the KVI [34|, whereas 225 Ra 
is best obtained from a 229 Th source. The radium ions 
from the cyclotron can be implanted in a solid AIF3 tar- 
get. RaF molecules can be produced from such a target 
by laser ablation. Similar molecules have been obtained 
in the rovibriational groundstate by the use of a cold 
buffer gas or supersonic expansion 35j, |36j . These pre- 
cooling methods greatly reduce the amount of photons 
that have to be scattered in order to reach the doppler 
temperature, which is a large advantage if the cooling 
cycle is not completely closed. 

A supersonic beam can be combined with a Stark de- 
celerator to bring the molecules to rest in the lab frame. 
A fraction of about 10~ 2 of the molecules in the su- 
personic expansion can be trapped. Heavy diatomic 
molecules such as RaF have not been decelerated to 
standstill yet; but a recent demonstration of an improve- 
ment to the Stark-deceleration technique [37|j indicates 
how this could be done. The efficiency of capturing the 
ablated molecules in the supersonic beam is however only 



about 10 ° 36, 38], which is prohibitively low for RaF. 



If the ablation is performed in a buffer gas cell, the 
cold molecules can be captured by a ma gnet ic trap [39j. 
The buffer gas has to be quickly removed [4Cf , and subse- 
quently lasercooling could be used to bring the trapped 
cloud to the lower temperatures that are required for a 
precision measurement. A promising number of 10 12 YbF 
molecules, produced in a single laser ablation shot, has 
been reported [3 51 ] but it remains to be seen how many of 
these can eventually be trapped and used for a precision 
measurement. 

To conclude, we have investigated the suitability of the 
radium fluoride (RaF) molecule to probe parity violation, 
which has so far not been detected in molecules. Explicit 
ab initio calculations of the parity violating nuclear spin 
dependent parameter W a and Franck-Condon factors for 
transitions between the electronic ground state and the 
lowest doublet excited state indicate that this open-shell 
molecule is a promising candidate for direct laser cooling 
and (subsequent) measurement of the nuclear anapole 
moment. We have also discussed possible approaches 



to creation, trapping and cooling of RaF molecules. In 
our opinion RaF is a valuable addition to the variety of 
molecules that may be used to probe parity violation, 
providing complimentary insight on the NSDPV effects. 
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FIG. 1: Morse potential fitting of the calculated potential 
energy points for the two energetically lowest-lying electronic 
states of RaF - the ground electronic state is 2 E 1/ / 2 and first 
electronically excited state is 2 il 1 / 2 - 
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